Conceptualization of seeded region growing by 
pixels aggregation. Part 1: the framework 

Vincent Tariel 
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Abstract — Adams and Bishop have proposed in 1994 a novel 
region growing algorithm called seeded region growing by pixels 
aggregation (SRGPA). This paper introduces a framework to 
implement an algorithm using SRGPA. This framework is built 
around two concepts: localization and organization of applied 
action. This conceptualization gives a quick implementation of 
algorithms, a direct translation between the mathematical idea 
and the numerical implementation, and an improvement of 
algorithms efficiency. 

Index Terms — Minkowski addition, seeded region growing by 
pixel aggregation. 



I. Introduction 

Many fields in computer science, stereo vision [5], math- 
ematical morphology [6], use algorithm which principle is 
Seeded Region Growing by Pixels Aggregation (SRGPA). 
This method consists in initializing each region with a seed, 
then processing pixels aggregation on regions, iterating this 
aggregation until getting a nilpotence [1][4]. The general 
purpose of this field is to define a metric divided into two 
distinct categories [2]: the region feature like the tint [1] 
and region boundary discontinuity [3]. The localization, the 
actualization and the organization of action in the interest zone 
is done manually. 

In this paper, we propose a conceptualization of localization 
and organization of action in the interest zone. First, a new 
mathematical object, Zone of Influence (ZI), is associated to 
each region to define a zone on the region outer boundary. 
Second as long as the algorithm is processed, the regions 
(de)grow, implying some ZI fluctuations. A procedure is 
defined to actualize efficiently the ZI such that some operations 
have to be computed only on the (de)growth neighbourhood. 
Third, in order to organize action on pixels belonging to the 
ZI, a system of queues (SQ) is defined to sort out each pixel 
belonging to the ZI depending on a metric and the entering 
time. Fourth, a library named Population has been created to 
implement these three first mathematical parts. An algorithm 
implementation using SRGPA with this library is quick and 
efficient. The translation between the mathematical idea and 
the numerical implementation is direct. Last, since this library 
has been optimized mathematically and numerically, each 
algorithm created with its is also optimized. In the article 
number 

II. Conceptualization of interest zone 

The basic idea is to define three objects: Zone of Influence 
(ZI), System of Queues (SQ) and Population. The algorithm 
implementation using SRGPA is focused on the utilisation of 
these three objects. An object ZI is associated to each region 



and localizes a zone on the outer boundary of its region. For 
example, a ZI can be the outer boundary region excluding 
all other regions. An algorithm using SRGPA is not global 
(no treatment for a block of pixels) but local (the iteration 
is applied pixel by pixel belonging to the ZI). To manage the 
pixel by pixel organisation, a SQ sorts out all pixels belonging 
to ZI depending on the metric and the entering time. It gives 
the possibility to select a pixel following a value of the metric 
and a condition of the entering time. The object population 
links all regions/ZI and permits the (de)growth of regions. 
A pseudo-library, named Population, implements these three 
objects. An algorithm can be implemented easier and faster 
with this library, fitted for SRGPA. The utility of this library 
will be developped in the article number 3 of this serie. 
This section is decomposed in four subsections: localization, 
actualization, organization, implementation. 



A. Localization 

The aim of the object Zone of Influence (ZI) is to localize 
the zone of interest where the algorithm is processing oper- 
ations. Its mathematical definition respects three conditions: 
simple, generic and numerically translatable. 
To understand the definition of this object, an example is 
given. Let assumed that it exists two types of lichens with 
two different colours, blue and red. At time t = 0, there is 
only a seed of each lichen on a desert island surrounded by 
water (see figure 1). At each step of time, each lichen grows 
on each pixel x belonging to its boundary if x doesn't belong 
to the other lichen or water. The region water never grows. 
The object ZI defines the localization of the growth for each 
region. 
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Fig. 1. Left figure: the white pixels are the desert island, the black pixels 
are the water and the red and blue pixels are the seeds of the two lichens. 
Right figure: after an infinite time, the lowest grey-level represents the blue 
lichen region and the middle grey-level represents the red lichen region. 



representing the red and blue lichens and the water. We define 
the object, ZI, for each region, Z* ed , Z f blue , Z l water . 



First the boundary can be defined as: 
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The symbol © means the Minkowski addition: X © V = 
^xgxV(x). V is an application of E to P(E), for example a 
ball. 

We add to this first definition the property that each pixel of 
the object ZI doesn't belong to the water and the other lichen. 
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This last definition defines correctly the initial constraints. 
In a general case, let X 1 = {Xq, . . . , X^} a set of regions. 
So for each region, we define a neighbourhood function, 
Vi = V or 1 and a restricted set N{ C N n such that each 
ZI associated to a region is defined as: 

Z\ = (Xt®Vi)\(\J X]) 

In the next papers[7][8], we will define differently the Vi 
and Ni to localize correctly the different zones of interest for 
each algorithm. This object ZI is generic and its mathematical 
definition is simple. The last property, numerically translatable, 
is presented in the next subsection. 

B. Actualization 

Firstly, we introduce two operators in the set theory, addition 
and subtraction. 

The addition, +, between two sets A and B, is the disjoint 
union: A + B = A l+J B 

The subtraction 2 , — , between two sets A and B, is: A — B = 

(A c + B) c 

We have defined these operators because they are useful in the 
numerical conversion of mathematical expressions. 
The SRGPA implies some fluctuations (growth or degrowth) 
of the regions. As the ZI depend on regions, the ZI have to 
be actualized. For the actualization, the simplest way is to 
calculate the fluctuation for each ZI. At each region pixel 
(de)growth, the numerical cost is high: the complexity is 
6(n) where n is the Zl/regions number. Numerically, the 
actualization has to be fast. 

A solution is to reformulate the link between ZI and the 
regions: the ZI at time t+1 depends on ZI at time t plus or 
minus a set defined in the neighbourhood of the fluctuation. 
This reformulation gives an actualization with a constant 
complexity 3 . 

! means an application of E to 

2 The addition is commutative, not the subtraction 

3 This section is quite technical and it is not necessary for this article 
understanding 



In this subsection, we suppose that we have one and only one 
region fluctuation A u between two steps of time: 



(yt e N,3\i e N n ) 



X] 



£+1 



Xj ± A 1 



VjeN n \i,x* +1 = xt) 



We see two parts in the expression of ZI: Z\ = 
myself other 

' * > TT* — * 

(X? m 0Vi)\( (J Xl Q ). The "myself part" depends 

only on the associated region. The "other part" depends 
on all regions. This is the reason why we decompose 
the actualization in two steps. First, we actualize the 



"myself part": X i 



£+1/2 



XI m ± A t then the "others part": 



= x !T /2 ± At ( see fi § ure 2 )- A11 the P roofs are in the 

appendice I. 

1 ) myself part: We suppose that there is only the fluctuation 
of X\ m between time t and t + 1/2. 

For all j different of i, Zj +1/2 is equal to Zj because they do 
not depend on X\ m . 

If it is a growth X^ /2 = X( m + A' then 

Z\ +1/2 = Zj + (A t ®V i )\Z t i \(\J XI) 

If it is a degrowth X*+ 1/2 =X t im - A 1 then 
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2) other part: We suppose that there is only the fluctuation 
of Xj a between time t + 1/2 and t + 1 
If it is a growth = X*+ 1/2 + A 1 then 



Vj G N n 



A* \ (Z* +1/2 ) 



if (i G Nj) A {Vj ^ 0) 



z t+i = z t+i/2 dse 



3 3 

The symbol A means the and in the symbolic logic. 
If it is a degrowth X^ 1 = X*+ 1/2 - A 1 then 

7 £+i/2, 
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3) Numerically: If we have a growth A 1 of the region X\, 
the ZI are actualized in two steps. 
First, the "myself part" actualization: Z^ 1 ^ 2 = Z\ 

V)\zt\( U xlo) 

. \/x e A* © Vi do 
if x £ Z\ then 

%xt U ^- then 
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■ x: 



4 Except the initialisation with seeds, each region fluctuation is a single 
pixel A t = {x} 
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Fig. 2. Two regions: X l r and Xj, two ZI Z l r = (X* V r ) \ (Ijulj) and 
= (X£ © y b ) \ (X* UlJ). Both series show the ZI actualization after a 
growth for the first serie and degrowth for the second serie. The first figure 
is the initial state, the second figure is the (de)growth of the blue region, the 
third figure is the actualization of the "myself part" of the ZI, the last figure 
is the actualization of the "other part" of the ZI. 



end if 
end if 
end for 

Second, the "other part" actualization: 
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Vj : x £ Z^ 1/2 do 
if i £ Nj A Vj + 

end if 



then 

- x\ 



end for 
end for 

If we have a degrowth A f of the region Xj, the ZI are 
actualized in two steps. 
First, the "myself part" actualization: Z\ 

^)\(4 + >^))\^) C 
. Mx e A 1 © Vi do 
if x £ Z\ then 

iixev^nx'tn 1 
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then 
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end if 
end if 
end for 

Second, the "other part" actualization: 5 
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Vx G A f do 

. Vj : (x e V" 1 ) n x^+V V do 

if ^ U ^ o 1 then 



if i i Z* +1/2 then 

5 The assumption that for all j, Vj is equal to or V , is crucial here. 



if i £ Nj A ^ ^ then 



end if 
end if 
end if 
end for 
end for 

The complexity of this actualization is time constant 6(1) if 
two data structures are implemented: 

L x (x) = {Wi:xeXi} 
L z (x) = : x £ Zi] 

At this step, the ZI localize correctly the zone of interest and 
there is a procedure to actualize the ZI efficiently. The next 
subsection will organize the action. 

C. Organization 

The general purpose of seed region growing is to define a 
metric [2] applied to the pixels belonging to the ZI. The priority 
order of pixel by pixel process depends on this metric and the 
entering time. At each iteration of pixel by pixel, a pixel is 
chosen because it respects some constraints on the metric and 
the entering time. The System of Queue (SQ) manages this 
organisation. 

The set B l is the couples to add to the SQ between the t + 1 
and t: 

B l = {(x,i):x^Z\Nx£Z t + x } 

The set C f is the couples to substract to the SQ between the 
t + 1 and t: 

C l = {{x,i) \ x £ Z\ Nx ^ Z l + X } 

The SQ is like a "factory store". There are: 

• the products, B l : , to store and to give a label, 

• the products, C £ , to remove from the store, 

• a customer that chooses product depending on its label 
and its entering time. 

To get the best efficiency in the "factory store" organisation, 
the three following points must be respected: 

1) Data structure: A data structure is a way to store data 
in a computer. A data structure depends on the operations to 
be performed, using as few resources both execution time and 
memory space, as possible. Generally [7], the data structure 
used is composed by n queues, Q l n = {ql, . . . , q^}, because 
the strategy is: 

• to store each couple (x, i) whose metric is the same in 
the queue (a queue is associated to a quantification of the 
metric and n is the metric range), 

• to organize in First In First Out (FIFO) "what comes in 
first is handled first, what comes in next waits until the 
first is finished". 

But the organisation can be a random access to handle ran- 
domly a couple [8]. 

2) The ordering attribute function: The aim of the ordering 
attribute function, S, is to assign each couple (x,i) of B l in 
the appropriate queue. The label k of the queue is given by 
k = 5(x, i) (see figure 3). NB: if the label k is equal to OUT 
then the couple is not affected to the set of queues. 



4 



S l (k, blue) = 5 
6'if, blue) = 3 



Fig. 3. Evolution of after the addition of the set B l 



N J ] Exlraction of ™ demeM ' here GlB 

1—2) Test if c 6 
Yes 

Fig. 4. How to act? l)Choose a queue 2)Extract a couple (x, i) 3)Test if 
x ^ Z\ 4) Yes, action, 5) No, come back to the first step 



3) The extraction function: Between time t and t + 1, 
all couple belonging to C l have to be removed in the set 
of queue. The numerical cost of this strategy is important. 
Another solution is to do nothing except when an couple (x, i) 
is extracted to a queue j at time t, then the test x G Z\ is 
operated 6 . If the test returns true, operations are processed on 
this couple. If the test returns false, then this couple is deleted. 
So the extraction function has two components: choose of file 
and test (figure 4). The general definition of the two objects, 
ZI and SQ, is finished and the next subsection introduces the 
algorithm implementation. 

D. The algorithm implementation 

In this section, a pseudo-library, called Population, is pre- 
sented 7 . This pseudo-library is composed with three objects. 
Each object is used in the elaboration of an algorithm. These 
three objects are: 

• the SQ that manages the "store factory". 

• the tribe that collects a restricted set and a neighbourhood 
to define after a ZI associated to a region. 

• the population that collects the set of regions/ZI and gives 
method for regions (de)growing. 

The implementation of this library is done using the C++ lan- 
guage. The conceptualization and the generic programming 8 
permit: 

1) an implementation of algorithms using SRGPA with less 
than fourty lines of codes, 

2) the application of these algorithms whatever the dimen- 
sion of the image (principally 2D, 3D) and the type of 
pixel/voxel, 

3) the optimization of all algorithms using SRGPA. Since 
the library has been optimized, all algorithms using this 
library will benefit from the optimization. 

, t—t' 

"When a couple (x, i) enters in the SQ, x belongs to Z i but is it still 

the case, now? 

7 Free software available at http://pmc.polytechnique.edu/~vta/population.zip 
8 An algorithm is not attached to a type of data but to some data properties. 



Object: SQ 

The aim of the SQ is to manage the "store factory". Its 
implementation implies the definition of: 

1) the ordering attribute function, 5(x,i), to give the label 
of the queue where a couple will be stored, for example: 

• 5(x, i) = 

• 5(x,i) = T(x) with X an image (an application of 

E to N n ) 

• ... 

2) the data structure of the elementary queue, (queue, 
linked list,. . .) 

3) the number of queues. 
This class provides three methods: 

1) select_queue(int i): selects the queue i (the customer 
selects a product with the label i). 

2) pop(): returns couple (x,i) from the selected queue (the 
customer gets a product with the label z). 

3) empty (): returns true if the selected queue is not empty, 
false else (the customer asks if there is still product with 
the label i). 

Object: Tribe 

The aim of the object tribe is to collect a restricted set 9 , N, 
and a neighbourhood, V to define a ZI associated to a region. 
This object doesn't provide method but is used by the object 
Population. 
Object: Population 

The object Population allows the creation of the set regions/ZI: 

{X\Z% = {(^,7Vo,y ,^),...,(^,iV n ,K,^)} such 
that Z\ = (X\ Vi) \ ( |J Xf). The initial state of the object 

Jem 

population is (A' t=0 , Z t=0 ) n = 0. So its methods permit the 
addition of region/ZI and the addition/subtraction of pixel(s) 
to a region. 

1) growth_tribe( Tribe tr=(V,N)): 

• creates an empty region/ZI, (X, Z) 1 with Z l = 

(^©y)\(U 

jeN 

• pushes this region/ZI at the end of the 
set of regions/ZI, Z t+1 ) n +i 

{(^z*) n ,{(x* +1 = 0,v;j\r,^ +1 )}} 

• returns the label of the created region/ZI in the set 
of regions/ZI. 

2) growth(Pixel x, int label _region/ZI): adds the pixel x to 
the region i, actualizes all the ZI and stores the set B l 
in the SQ using the ordering attribute function. 

3) growth(Set A, int label_region/ZI): adds the set A to the 
region i, actualizes all the ZI and stores the set B l in 
the SQ using the ordering attribute function. 

4) degrowth(Pixel x, int label_region/ZI): subtracts the pixel 
x for all regions i if x belongs to Xi, actualizes all the 

Numerically, the object restricted set, N, has two fields: 

• cl-w (all or without), a boolean 

• L, a list of integer 

if 

_ \ N n \ L if a.w = true 
| L else 



ZI and stores the set B l in the SQ using the ordering 
attribute function. 

5) Z()[x] gives the list Z x = {Vi : x £ 

6) X()[x] gives the list X x = {Vi : x G 
Algorithmic in general 

An algorithm is usually presented with: 

1) Definition of SQ: the ordering attribute function 5, the 
elementary queue q and the number of queues 

2) Definition of the different tribes 

3) Creation of regions/ZI using these tribes 

4) Initialization of the regions/ZI with seeds 

5) The growing process: 

Select a queue 

while the selected queue is not empty do 

Extract a couple (x,i) from the selected queue 
Process some operations like 

(de)growth_pixel(x,i) 

end while 

6) Return result 

The algorithm 1 (see figure 1 and 5) is an example of these 
steps linking. 



Algorithm 1 the growth of lichens in a desertic island 
(geodesic dilatation) 

Require: I, S , V //The binary image (I(x) ^ means x 
belongs to the desertic island and else belongs to the 
water), the seeds, the neighbourhood 
II initialization 

System_Queue s_q( 5(x, i) = if I(x) + 0, OUT else, 
FIFO, 1); //A single FIFO queue such as the pixel be- 
longing to the water are not pushed in the SQ. 
Population p (s_q); //create the object Population 
Restreint TV = N; 
Tribe actif(V, N); 
.\fsi e S do 

int reLtr = p.growth_tribe(actif); //create a region/ZI, 
(XI Z\) such as Z\ = {X\ © V) \ ( U x j) 

p.growthfe, reLtr ); 
end for 

// the growing process 

s_q.select_queue(0); //Select the single FIFO queue. 
while s_q.empty()==false do 

0, i) = s_g.pop(); 

p.growth(x, i ); 
end while 
return p.XQ; 



III. Conclusion 

In this paper, we have conceptualized the localization and 
the organization of seed region growing method by pixels 
aggregation. 

In the conceptualization part, we define two objects and one 
procedure to make possible the creation of the library, called 
Population. The first object, zone of influence, is associated to 
each region to localize a zone on the outer boundary region. 
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Fig. 5. The first figure is the desertic island, the second figure 
is the set of seeds, the next figures are the geodesic dilatation such 
as each colour represent a specefic region. A video is available at 
http ://pmc .poly technique . edu/ ~ vta/geodesic .mpeg 



The second object, the system of queue, organizes the queue of 
pixel by pixel aggregation around a concept of "factory store". 
The procedure is defined to actualize at time constant the zones 
of influence after a fluctuation of a region. This procedure and 
these objects permit to this library to be numerically efficient 
and to implement algorithm faster. 

This first paper will be follow by others. In the paper part 
2, we will present a procedure to define a boundary between 
the region such as its localization does not depend on the 
seeded region initialisation order. In the paper part 3, this 
library will be used to implement a wide range of algorithms. 
In the paper part 4, a simple, generic and robust method 
will be proposed to extract the components from experimental 
tridimensionnal images of granular materials and porous media 
obtained by X-ray tomography using algorithms coming from 
SRGPA. In the paper part 5, the cutting of the porous media 
in elementary pore will be done following two different 
conventions: morphological and topological. In the paper part 
6, we will propose an efficient procedure to reconstruct a three- 
dimensional (3D) medium using morphological information: 
chord length distribution and two-point correlation. The effi- 
cient is due to the localization of the interchange procedure in 
the boundary of the phases. 

Appendix I 
Proof of the acualization 

A. Growth of the myself region 

If there is only the growth of the myself region: X 1 ^ 2 = 
Km + At ^en 

^• +1/2 = z? + (A*e^)\^\(U x lo) (D 

j£Ni 

Proof: 

zt +1/2 = (4r /2 ©^)\(U4: l/2 ) 

jeNi 

As we have only the growth of the myself region: X^' 2 = 
Km + A* then 

Zl +1/2 = ((Xl )m + ^)©T/,)\(U X& 
As (AUB)®V = (A©V)U(5©V) so withi?! = (J Xj 

Zj +1/2 = ((Xl m (BV i )U(A t (BV i ))\Rl 
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We know that A U B = A + (B \ A): 



yt+l/2 



{{x{ m © Vi) + (A* © Vt \ (xl m © v))) \ R\ 



and (A + B) \ C = (A \ C) + (B \ C), then 



7 t+l/2 _ 



(Xt © V) \ R\ + (A* © V \ (X* m © Vt)) \ R 



Wchayc Zj = (Xl m ®Vi)\Rt 

z t+i/2 = Z t + {A t (BVi) \ {X t m(BVi) \ R t i 

As (A\B)\C = (A\(B\C))\C then with Z\ = (Xf , 
V)\Rt 



7 t+l/2 



Z\ +1/2 + {A 1 © V;) \ Z\ \ Ri 



B. Degrowth of the myself region 

If there is only the degrowth of the myself region: X\ 
Vm +^ then 



^• +1/2 = Zl-{{A^V % )\{X^ /2 ®V % ))\{Zl)\2) 
Lemme : 

A © V = ((A U B) © V) \ ((B © V) \ (A © V)) 
First, demonstrate that A = (A U B) \ (B \ A) 



A = 


A\J% 


A = 


iu(Bn B c ) 


A = 


(iuB)n(iu B c ) 


A = 


(A U B) n (B c U A) 


A = 


(A U B) \ (B c U Af 


A = 


(AuB)\(BnA c ) 


A = 


(AuB)\(B\A) 



Switching A by A © V and 5 by B © V in this last formul, 
we get: 

A © V = ((A © V) U (B © V)) \ ((5 © V) \ (A © V)) 
As (A © V) U (B © V) = (A U 5) © V, we find 
A © V = ((A U B) © V) \ ((5 © V) \ (A © V)) 
Proof: 



Z 



£+1/2 



«r /2 ®^)\(U4: l/2 ) 



As we have only the degrowth of the myself region, X 1 ^ 2 = 
Xl m - A* with R\ = U Xl„ then 

z* +1/2 = ((xt^-A^e^XRi 

By the lemme, switching A by (Xj m — A 1 ) and i? by A*, we 
have 



zt+l/2 



= (A 1 \A 2 )\Rj 
with Ai = ((X* m - A*) U A 1 ) © ^ 
with A 2 = (A t ®V i )\ {{X\ m - A 1 ) © Vi) 



Like ((X* m -A*)UA* = X* m so we obtain ^ = Xf, m ©Vi 
As X* m - A* = X<+ 1/2 so A 2 = (A* © V) \ © V,). 

We have 

^• +1/2 = {A 1 \A 2 )\R\ 
with Ai = X? m © Vi 
with A 2 = {A 1 © Vi) \ (X*+ 1/2 © V) 

As (A\B)\C = (A\C)\B, substituting A by Ai, B by 
A 2 and C by i?- in this last formula. As A 1 \ R\ = we 
have 

Z* +1/2 = Z|\A 2 

with A 2 = (A t (BV i )\ (X*+ 1/2 © Vi) 



t+i/2 _ As A\B = A - (B\ A c ), we have: 

z t + i/2 = z t_ {{A t &Vi) \ {x ^/^ Vi))X{z ty 



C. Growth of the other region 
If there is only the growth of the other region: X\\ 

x t+l/2 + At ^ 



-A* \ (Z* +1/2 ) c if i 6 JVj A V} 7^ (3) 



= ^ 



£+1/2 



else 



Proof: 



Let = 0, then Vt : Zj = 0, in particulary Z] +1 = Zj +1/2 . 
Let assume now that Vj = V. 



Z^= (X^(BV)\U keNj Xt+; 
By the commutativity of union, we have: 

Z J +1 = ( X £m © V) \ {U keNAi X\+ l U X* +1 ) if i G ^ 



Z t+1 = 



else 



We suppose i G iVj. We have only the growth of the other 
region: X™ = X^ 2 + A* 

Z] +1 = (X^ 2 © V) \ (U keNjV Xt+ 1/2 U (Xt +1/2 + A f )) 
A + B = AD B and by the commutativity of union, we have: 



(C /2 ®^)\( u ^C 2uA ') 



As A \ (B U C) = A \ B \ C and Z* +1/2 = (Xj+ 1/2 © V) \ 
(Ufc eJVj X^ 1/2 ), we have 



zr 



As A \ B = A - B \ A c , thus 

Z t + 1 = Z t+l/2 _ A t ^ {Z t+l/2y 



7 



D. Degrowth of the other region 

If there is only the growth of the other regiomX^ 1 = 



X 



£+1/2 



A 1 then 



Vj G N n ^ 



zj +i 



7 t+l/2 



(((A*\(xi + ] /2 e^) c )\( U x*+ 1 ))\zj +1/2 ) 



if (< € Nj) A (V^ 0) 



z ,+i 



Zj +1/2 else 



t+i 

i 



Proof: 

Let Vj = 0, then V£ : = 0, in particulary Z 
Let assume now that Vj = V. Let j in N n 

z* +1 = (x£>v)\u fceA ,,x«+ 1 

By the commutativity of union, we have: 

= © V) \ (U keNAt X^ U X^) 

r?t + l _ gt+l/2 



rjrt+l/2 



[1] 

[2] 
[3] 

[4] 
[5] 

[6] 
[7] 
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si i G Nj 
ifi <£ Nj 



We suppose i G Nj. We have only the growth of the other 



region 



:X^=X 



t+l/2 



A* 



= «+ m 1/2 © V) \ (U keNjV X^ 2 U (X*: 1 / 2 - A*)) 



3 

As A 

Z^ 1 



B = A fl B c , we have 



(4+ 1/2 e v) \ (u kGNAi xt + 1/2 u (x<+ 1/2 n (^) c )) 



As A U (B n C) = {A U B) n (A U C), thus 

© V) 



(v t+l/2 



\(u keNj xl!o 1/2 n (u keNAl xtt 1/2 u W)) 



As A \ (B n C) = (A \ B) U (A \ C), thus: 



zr = 



Z t j +1/2 UA 1 



(X-^ /2 ®v)\(u keNA ^ k 



with A\ 

As A \ (B U C) = (B c DA)\C, we have 



Xlt 1/2 U (A*) c ) 



Ax = (A* n (x*t! /2 © V)) \ (u keNjV xl^) 

A, = (A* \ (X*+ 1/2 © Vf) \ (U k€NAl Xi+ 1/2 ) 



t+l/2s 



A 1 = (A t \U, eiVAi X*+ 1/2 )\((X 



.t+1/2 
j,m 



®V) C ) 



X\ + 1/2 - A\ so A 1 n X*+ 1/2 = and A 1 



As Xl Q 

A* \ XI+ 1 . As (A \ B) \ C = A \ (B U C), we have 



t+1/2 ux^))\((x^ /2 ®vT) 



Ai= (A* \ (U keNAi X kt0 

M = (A* \ (u keNj xj^) \ ((xj+ 1/2 © vy) 

m = ( At \ (4- 1/2 © y) c \ (Ufc^ft 1 ) 



The last step is: zj +1/2 U A 1 



Z* +1/2 + A!\Z* +1/2 . 
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